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We show that incorporation of strong electron correlations into the Kohn-Sham scheme of band 
structure calculations leads to a modification of the standard model of the lanthanides and that this 
procedure removes the existing discrepancy between theory and experiment concerning the ground 
state properties. Within the picture suggested, part of the upper Hubbard /-band is occupied due to 
conduction band- /-mixing interaction (that is renormalized due to correlations) and this contributes 
to the cohesive energy of the crystal. The lower Hubbard band has zero width and describes fermionic 
excitations in the shell of localized /-s. Fully self-consistent calculations (with respect to both charge 
density and many-electron population numbers of the /-shell) of the equilibrium volume Vb and the 
bulk modulus of selected lanthanides have been performed and a good agreement is obtained. 
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For strongly correlated systems, state-of-the-art band structure methods fail to give an accurate description. The 
reason for the failure is well-known: the strong intra-shell Hubbard repulsion, U , is underestimated in the calculations. 
If the /-/ interaction is sufficiently strong, the /-shell forms a localized multiplet (Russel-Saunders coupling) and can 
accurately be treated as atomic like. This is the basic idea behind the so-called standard model of the lanthanides and 
many of the physical properties seem consistent with it, for example the Curie^Weiss behavior of the susceptibilityfil, 
the gross features of the equilibrium volumesa, and the structural propertieg3. When using ah initio calculations 
employing the local density approximation (LDA) and treating the /-electrons as core states, a too large equilibrium 
volume, Vo, is obtained (for all other elements, LDA underestimates the equilibrium volume). On the other hand, 
when the 4/'s are treated as itinerant electrons, Vq is much too low compared to experiment and the observed localized 
moment (Curie- Weiss law) is absent. Calculations by Delin et alu show that treating the 4/-electrons as core states is 
a very good approximation for the late lanthanides, but for the lighter lanthanides the disagreement between theory 
and experiment is gradually increasing when one goes from heavy to the lighter lanthanides. Thus, even though 
the standard model is essentially correct, it has to be modified slightly to better describe the cohesive properties of 
the lanthanides. Such a modification requires a method combining an ah initio band structure calculations (AIBSC) 
and a roany-body approach. A variety of such methods are developeid. The ideas of orbital polarizationQ, LDA-I-U 
methodtl, Hartree-|Ebck type of approximations, method of phase shiftu and more advanced methods, like three-body 
Faddeev equationsEHj, have in some cases been successfully appMed to systems with strong electron correlations (SEGX 
New approaches: on the basis of the Gutzwiller wave functionB, dynamical mean fieldlia and the LDA+ -I— approacnllU 
have also been used for the description of different SEC systems. In many cases the AIBSC-based density of electron 
states is used as input to a model calculation using a non-self-consistent fixed lattice parameter and Hubbard U . 
However, since a large value of Hubbard repulsion restricts the available electron phase space, wave functions and, 
correspondingly, the rriartrix elements of mixing and hopping, changes. Recently an attempt to describe SEC by means 
of slave-boson methodllHl (i.e. Hubbard U — oo) build into the_scheme of linear muffin-tin orbital method (LMTO) 
was made. The author found that the results are unsatisfactoryll3 and had to use a scheme different from the standard 
slave-boson self-consistency schemes. However, in all these attempts the infiuence of SEC on the cohesive properties 
has not been questioned yet. Our approach to improve on the standard model is based on the following physical 
picture: the Hubbard repulsion is sufficiently large to form large energy gaps between the (n — 1)-, n- and (n -I- 1)- 
electron configurations not only in the heavy but also in the light lanthanides . Near the atomic limit all /-electrons 
are in the strong coupling regime and hopping and mixing interactions lead to the formation of a lower Hubbard sub- 
bands (LHB, |n — 1) — > |n) transitions) and a upper sub-band (UHB, \n) \n+l) transitions). In our case, however, 
the lower transitions are substantially below the bottom of the conduction bands (Ce is an exception) and, hence, they 
do not experience any mixing interaction. Hubbard U increases with atomic number, due to the more localized nature 
of the 4/'s, thus the LHB goes down in energy while the UHB goes up and therefore it experiences smaller mixing and 
is populated with fewer states. These many-electron local excitations correspond to localized core electrons within 
the one-electron picture. The upper transitions are in the energy region where the delocalized electrons do exist. 
Therefore, the overlap and mixing interaction may delocalize these transitions, forming an UHB, filled only by a small 
amount, 77, of electrons. This picture corresponds to the multi-orbital Hubbard- Anderson model. The total occupation 
number can be written as Uf = n + rj, where n is an integer, which is determined by the nuclear charge valence, and 
< 77 < 1. Thus, most of the /-occupation is in the LHB. This is consistent with the observed Curie- Weiss behavior 
of the high-temperature (above Tc) susceptibility. In the lanthanides the width of the appropriate /-sub band is much 
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smaller than the Hubbard band gap. This allows one to formulate an-efFective Schrodinger (or Dirac) equation (SE) 
for two different energy regions. As is known from Hubbard's modelstS, the single-electron spectral weight at J7 = 
is split by the Hubbard repulsion [/ to a set of weights, corresponding to an electron transition from one ion state 
to another. These weights renormalize the matrix elements of the hopping and mixing interaction, which in turn, 
cause a narrowing of the bands due to correlations. The theory is developed for a non-orthogonal basip_set, as for 
instance used in the linear muffin-tin orbital method in the atomic sphere approximation (LMTO-ASA)Efl. Here we 
apply the formalism only in the lowest approximation which requires less complicated changes for methods used in self 
consistent ab initio calculation. Compared to conventional considerations of the periodical Anderson model (PAM) 
which do not include structure-dependent vertex corrections we solve the Schrodinger equation for the electrons in 
fields generated by the nuclei and electrons in every iteration, with the charge density obtained from a system of 
equations for Greens functions (GF) and spectral weights, thus renormalizing self-consistently the matrix elements of 
the PAM. On the other hand compared to the normal LDA approach wc solve self-consistently additional system of 
equations for many-electron spectral weights and obtain a bandwidth which is reduced by correlations. 
Technically our method, briefly described, consists of the following steps: 

1. The total many-electron Hamiltonian is written in terms of an LMTO basis set. 

2. The /-electron operators are transformed into a Hubbard representation. 

3. Analytical calculations are performed for the-electron Greens functions for the corresponding PAM in the LMTO 
representation in the Hubbard-I approximationll3 (strong-coupling regime). 

4. Renormalizing factors to the f-{spd)-h\ocks of the LMTO Hamiltonian and overlap matrices are introduced which 
are derived from a comparison of the frequency and overlap matrixes. They arise due to strong correlations in the 
system for Greens functions. 

5. Self-consistent, ab initio calculations of the system with SEC ajCe performed. 

In practice, the main effect of these steps consists of the separationlla of the /-electron system into two subsystems and 
the correlation driven narrowing of the /-bands in the conventional LMTO Hamiltonian. This, as will be demonstrated 
below, leads to an improved description of the bonding of the lanthanide metals. In order to describe our theory in 
more detail, the secondary quantized full many-electron Hamiltonian, Ti., is written in terms of an LMTO-basis set. 
An unperturbed Hamiltonian and a perturbation is formulated as Ti. — {T + Vne + Vlda) + {Hee — Vlda) = Tio + W- 
The main effect of the /-part single site correction from Coulomb interaction, 

'^U 2 ^ ^ ^mim2m3m4,fmiijfm-,a-'f'm3C7'frn4,cr ^ '^^^ ^1=3. a ^ ] fm,crfmicn (1) 

is to localize part of the /-spectra, similar to a so called LDA+U solution. Let us denote all non- /-operators in the 
LMTO representation as CkL- Then the LMTO part of the Hamiltonian, J2 -ffLL'(k)aj^^akL', can be written as a sum 
of s,p,d, and /-electrons, and their mixing interaction, 

nmia kL,L' 

tnm '/nmicr/mmjo- + [^^L,r?i,(T (k)e " C^i^^fnmicr + h.c], (2) 

nm,m;crmj k,n, L,mia" 

where inm"' = X^k [-^™!<^,™'o-(k) — <5mi,m'e/] is the hopping matrix element between /-orbitals and e'^ = 

J2k^-mi<T,mi<T0<)- Since the LMTO-ASA basis is non-orthogonal with an overlap matrix (xnilXmi') = OnL.mL', 
the anticommutator is {anL,a|nL'} ~ (C^^)nL,mL'- Therefore, the equation for the bare fermion Greens function, 
G„L,mL'('^) = -{T[anL{t)al^j^,{0)])^, for this Hamiltonian is 

[ujOnL,miLi — ^^nL^miLi ] Gmi ,Li ,mL' (t^) — ^n,mSL,L'- (3) 

Thus, we find (by construction) in the square brackets the conventional secular problem for the LMTO-ASA method 
(or any LDA based method). Introducing a full set of states in an orbital representation jFo) = |0), jFi) = I7), IF2) ~ 
ITjT'))''' (7 denoting a single particle /-state, mjcr) allows to express any single-site operator A in terms of X- 
operators, A = J2r r'i^l^l^')-^^'^ ' ^^"^ — \p){l\ where {X^'^) = Nr is the occupation number for the state 
F. Now we can express the zero /-Greens functions, for the LMTO problem, in terms of X-Greens functions as 
G^"\uj) = \if'r)a\^PaDa{u;), whcrc Daico) = and an X-operator GF G^'^ = P^Da. Here a = (F„,F„+i) 

(a = (F„+i,F„)) are possible Hubbard transitions, Pr„,r„+i = A^r,, + A'^r„+i, Ar„+i,r„ ~ ^r„+i — Er„ — fJ., -Er„ — 
ne° -|- Un{n — l)/2 and /i is the chemical potential. 

The /-electron operators can be expanded in the Hamiltonian Ti. into X-operators for the set of many-electron 
states described above. Then we arrive at an Anderson like model 
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n^Y. ^^^n^ + E ^L.L' (k)4^CkL' + (4) 
nr k,L,L' 

k,n,/i,a 



Here — iy(/n^i)°(/mi/)^*nm- I^oi' the derivation of the equations for the Greens functions, CJ^'^^^ = 

-{T[cnL{t)X^^^{0)]) and g'-^^'^ = -{T[X™ {t)X^^^ (0)]) , (which will be compared to the equation for the regu- 
lar Kohn-Sham (LMTO) problem) we need the anticommutators {c,X}, we find {cnL, X^,^} = 

X^m^riTa ^nl nitj.iff^)T"i,r2{Xni^^ J X^,'^}. Physically, this implies that the /-part of the conduction-band wavefunc- 
tion, coming from other sites, experiences strong correlations on the reference ion. The Hubbard-I as well as mean 
field approximation correspond to {X^^'^^ , X^ ~ -Pr',r^r',r2'^r,ri • Introducing renormalized "Fermi" operators 
= X^''^ / y/Pr'^r- We see that the matrix form of the GF, Q'^^^\ with i,j = c,X, coincides with the GF of the 
Kohn-Sham equation, G, with the only difference being A^.^ = J2r' r ^L,fj.iffi)r',r\/Pr',r, where A = H ot 

O. Thus the most essential, technical part of our theory leads to an effective secular equation which, when diagonalized 
gives eigenvalues with many-body corrections. 

As an example let us consider praseodymium metal, which has two localized /-electrons. Therefore, |ej| has 
a value which gives an energy minimum for n ^ 2. For briefness we will use the following notation for the /- 
orbitals: j = {rrii = 3,i) = 1, (m; = 2, i) ee 2 and 7 = for 7 7^ 1,2. r(i2) = (12), r(i2i,) = (Uv). Let us 
now consider a simple polarized solution, where the /-orbitals with 7 = 1,2 are fully occupied, while the rest of 
the /-electrons occupy the /"^-states, jTs) = |l,2,z/). When mixing is absent, n/ — 2 and ^''"^2) — 1 whereas all 
other population numbers, N^, are zero. Then = £^2(12)1^^7=2 = Di(^i2) and Gj=i, = £^(i2)(i2i/)- We denote 

^(12)2 = ^(12)1 = ^1, and A(i2^)(i2) = A2. In the hmiting case, A2, A3 00, no /'s are present in the conduction 
band, and therefore this corresponds exactly to the limit of the standard model. The /^-bands are slightly above the 
Fermi energy, and mixing interaction transfers /-character into the s,p, d-electron states. In a general case it is, of 
course, impossible to write equations for Nr in terms of only orbital GF's, G-y, because equations for higher correlation 
functions are needed. However, if for simplicity we assume that these 12 !^-bands are occupied symmetrically, i.e., 
^^(12!.) = (1/12) Af(i2i.'), we immediately find 
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= 2 ■ Ar(i2) + 3 • ^(12.) = 2 -f 77, 



^(12)+E.^(12.) = 1, 



(5) 



i.e., N{i2v) = ^7/12 and A''(i2) = 1 — rj (for the other rare-earth elements a slight modification of these expressions are 
needed). Thus, renormalization of the population numbers leads to the form of /-locators (single-site, k-independent 
part of the GF) 

r r , ^ J r 1-11^12 , 1W12 

Qi=Q2 = T- + T- and = + — , (6) 

uj — Ai uj — A2 uj — A2 uj — A3 

where A3 = A(i2iyiy')(i2iy')- The center of the non-rcnormalizcd LDA-/-band is in a; = A2 and Q^^^ = Qp^^ — 
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Qu ^ l/C"^ ^ ^2)- Let us ignore for the moment the non-orthogonality. For a diagonal effective hopping tjCkjLu), 
of any origin, we have G-y — [(Gq"*^)^^ — tj]~^. Therefore, the expressions for the local, /-part, of the self-energies 
are = {Gil^^)~^ — Uj^^y'i^ ~ S^^~f>)n^r — {g'"q!^)~^. Since all poles in are well separated, one can formulate 
two effective Schrodinger equations near w ~ Ai and lo ^ A2. Renormalization factors (1 — {^)uj=Ai) are, of course, 
given by spectral weights (numerators) of the GF's in these poles. Since 77 is small the bandwidths of the upper 
Hubbard sub bands for 7 = 1,2 are also small and, being above the Fermi level, they are empty. Within the scenario 
considered the lower sub bands have weight slightly less than one, (1 — 77). For the bands v, A3 is far above the Fermi 
energy. Thus the upper sub bands are empty and the spectral weight in the pole lu = A2 is (1 - i3£i). Therefore, 
the effective Schrodinger equation for bands 1 and 2 at w ~ Ai = A2 — C/ gives uj — (A2 — C/) — (1 — J])ti(k,uj) = 0, 
whereas for the bands at w ^ A2 we find w — A2 — (1 — llr7/12)i^(k, lu) — 0. 

Since in the case of LDA calculations with delocalized /-electrons the center of the /-band is ^^^2 we conclude that 
the potential for the Ist-and 2nd orbitals is shifted down by U (similar to the LDA-f U - resultjl3 and self interaction 
correction (SIC)-theorjil3 but without factors n^). The next conclusion is that the described factors should be to 
introduced into the LMTO overlap and Hamiltonian matrices. Thus, we calculate the charge density and a self 
consistent potential, corresponding to the given j], where i] comes from the equation of self-consistency 

r; = -i Y /(-)«r'^^'')'^^^'^'^^H- + *^,k), (7) 
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since the parameter Pr'.r = 1 — renormalizes Ti. and O. Here g^^^'^^'^''''^^^'^'^^^ jg obtained from the transformed 

secular matrix for the normal LMTO-GF's, {OE — H)^^, via /. For Pr we find the following picture. Two of the 
f-orbitals are located at an energy -U lower than the remaining /-states and below the bottom of the conduction 
bands. Therefore, we treat those /-states as core states removing the somewhat arbitrary choice of a value for the 
Hubbard U. Compared to LDA-I-U and SIC methods, in our approach the spd-f hybridization is reduced by the 
renormalizing factor a/P. The theory is generalized for the other lanthanides straightforwardly. 

Before giving numerical examples of cohesive properties (e.g. equilibrium volume) obtained from our method we 
describe how the total energy was calculated. As is well-known, the equation for the exchange-correlation potential 
Vxc derived by ShanJl3 connects it with the exact self-energy of the electron system. Making use of a strong-coupling 
(SC) perturbation theory (PT) developed by us we have performed the analysis of the contributions to the self-energy 
and found that a there exist a one-to-one correspondence between the graphs for the self-energy in the standard 
weak-coupling (WC) PT (in e.g. Jiie random phase approximation) and the sequence of graphs in SCPT. Using 
the facts that: a) Sham's equatiorO should produce a Vxc corresponding to the patticular SE choice; b) established 
correspondence between the WC and SC PT graphs; c) the statement of KotaniE2l, that the static random phase 
approximation reproduces well the DFT-LDA calculation with the standard choice of v^c, we come to the conclusion 
that, within the approximation chosen, the original LDA functional form for v^c can be used. However, the charge 
density is calculated from renormalized fermions which differ from the standard fermions by the renormalizing factors 
•\/P. Hence the charge density,pp(r), also depends on t/P. The mixing and hopping in the effective Kohn-Sham 
equation should be renormalized by the same P as described above and a part of the renormalized /-electrons 
should be described as localized (lower transitions ) and the rest (coming from the upper transitions ) as delocalized. 
The separation is caused by self-interaction correction automatically generated in PT. Finally, the equation of self- 
consistency for parameter P (Eqn |^) should be added. 

We have hence calculated the total energy using the regular LDA total energy expression, using the electron density 
pp(r), i.e. ELDA{pp{r)^ where the parameter P which comes from electron correlations, is found self-consistently 
from additional equation and regulates distribution of spectral weight between lower and upper parts of energy. 
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FIG. 1. Fig. A shows the equihbrium volumes, in experiment, LMTO-ASA (using the standard model), corrected 
LMTO-ASA and shifted theory, for some selected lanthanides. The basis-set shifts used are (in a.u.) for Ce Ai?sms =-0.157, 
-0.185 for Pr, -0.184 for Nd, -0.176 for Pm, -0.172 for Sm and -0.003 for Tm. Fig. B shows the bulk modulus for the same 
elements. The slope of the lanthanide correction is reduced due to the contribution to the cohesive energy from the delocalized 
/'s. 

In Fig. we show the equilibrium volumes of some selected lanthanides, from experiment, uncorrected Kohn-Sham 
calculations (LMTO-ASA) and from our corrected theory. Note that the many-body corrections are presented together 
with shifted data due to an incomplete basis set used in the LMTO-ASA calculations. This shift was calculated by 
comparing (LDA) volumes obtained from a full-potential-(FP)-LMTO, multi-basis set calculation with (LDA) volumes 
from LMTO-ASA calculations, ARs^s = Rq^ - R^^'^'^~^^'^ ■ An overaU good agreement is obtained, and for the 
heavier lanthanides the renormalization correction becomes vanishingly small (as seen for Tm). The only parameter 
of our model,|-t|he position of the UHB with respect to the Fermi level, was taken from bremsstrahlung isochromat 
spectra (BIS)Eil. In principle this value can be found from a super cell calculation changing the /-occupancy at one 
site. For cerium, A„^„+i is ~3.5 eV above the Fermi level, and shifting this to 3.5± 1 eV, only affected the equilibrium 
volume with ±0.5%. For samarium the position of the UHB is ~0.6 eV, and changes in the position with ±0.1 eV, 
changed the volume with ±0.8%. For the other elements changes were of the same order. Therefore, we conclude 
that the results are rather insensitive with respect to the position of the UHB. Fig. |i|B shows the bulk modulus. Note 
that this value is quite sensitive to the fitting procedure (pressure-volume fit) and to the structure, thus, it should be 
considered only as an approximate value. l— , 

The physical picture obtained is consistent with previous theoretical work on Fermi surfacesEZI where a better 
agreement between theory and experiment was obtained when including, ad-hoc, the 4f states as valence electrons 
even though part of these 4f states were treated as occupied_core states. It is also consistent with results obtained from 
self interaction corrected density functionals on Pr metallj. Unlike the method of orbital polarization the present 
theory is a true many-body correction, and thus presents a different physical approach. The decisive parameter in our 
case (unlike to the LDA+U method) is not the absolute value of the Hubbard U but the position of the transition 
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E{Tn+i) — E{Tn) with respect to Fermi energy. We have taken this parameter from experimental datau and this led 
us to an improvemenjtin the calculated cohesive properties. First loop corrections, affecting the position of the levels, 
are small in our caseEa. Fluctuations coming from charge excitations are expected to alter the situation at non-zero 
temperature, however, this remains an open question. 

We have shown that many-body corrections to the standard model improves the description of the cohesive prop- 
erties of the lanthanides, and explains the observed discrepancy for the equilibrium volumes for the light lanthanides. 
The only input to the present theory is the position of the upper Hubbard transition, Ar+i,r, which was taken from 
bremsstrahlung isochromat spectra (BIS) data. We found that a moderate change in energy of this transition gives 
only a minor change in the equilibrium volume. The calculated volumes are in good agreement with the experimental 
ones, which means that the excited states contributes to the ground state properties. The presented theory is designed 
to deal with strongly correlated electron systems and is integrated in an electronic structure method, and made fully 
self-consistent. It describes an additional contribution to the chemical bonding from a fraction /-electrons hybridized 
with the conduction electrons. Compared to for instance the LDA-f U, orbital polarization or SIC methods the present 
theory involves a decrease of mixing due to strong correlations, causing the /-bandwidth to be reduced. 
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